Perfect transmission and highly asymmetric light localization in photonic multilayers 
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General principles for the existence of perfect transmission resonances in photonic multilayer structures are 
formulated in terms of light interference described by recurrent Airy formulas. Mirror symmetry in the multi- 
layer is shown to be a sufficient but not necessary condition for perfect transmission resonances. Asymmetric 
structures displaying perfect transmission in accordance with the proposed principles are demonstrated. A hy- 
brid Fabry-Perot/photonic-crystal structure of the type (BA)*(AB)*(AABB)"' is proposed, combining perfect 
transmission and highly asymmetric electric field localization. Strength and asymmetry of localization can be 
controlled independently, to be of use in tailoring non-reciprocal behavior of nonlinear all-optical diodes. 

PACS numbers: 42.25.Bs, 78.67.Pt, 42.25.Hz, 42.65.Pc. 



I. INTRODUCTION 

Probably the simplest case of inhomogeneous media, pho- 
tonic multilayers are a good testing ground for structures with 
complex geometrical properties such as aperiodic long-range 
order (see, e.g., 0JJ and references therein). Indeed, the avail- 
ability of simple, cheap, and reliable computational methods 
often makes it possible to relate geometrical and optical prop- 
erties in an explicit manner. To name a few examples, scaling 
and self-similar features in optical spectra of quasiperiodic Fi- 
bonacci J2-01 and fractal Cantor multilayers were re- 
cently found to result from geometrical self-similarities of the 
underlying structure. It is even possible to formulate general 
relations for spectral properties of structures with arbitrary 
layer arrangement [HI]. 

One of rather intriguing properties of aperiodic multilayers 
is the appearance of perfect transmission resonances (PTRs) 
in the optical spectra, i.e., frequencies for which the multilayer 
has transmittance exactly equal to unity (\T\ = 1). It is known 
that multilayers with mirror symmetry (e.g., Cantor) com- 
monly exhibit PTRs while those without it (e.g., Fibonacci) 
usually do not: transmission peaks in such multilayers, even 
if they look "perfect", really have \T\ < 1 (see Fig.Q]). Sev- 
eral accounts II 91 41211 report PTRs if a Fibonacci structure is 
symmetrized and show that perfect transmission is explicitly 
related to mirror symmetry j 1 1(1 . However, more recent results 
show perfect transmission in asymmetric multilayers based on 
periodic lfT3tl . Fibonacci Hill , and Thue-Morse IU5ll geometry. 
This suggests that mirror symmetry is sufficient but not nec- 
essary for PTRs. 

Such PTRs in asymmetric structures are promising in de- 
signing non-reciprocal optical devices such as nonlinear all- 
optical diodes llol . Indeed, an associated spatially asymmet- 
ric light localization at resonance (see IU3l[l5lO induces a non- 
reciprocal nonlinear optical response, while perfect transmis- 
sion assures that reflection losses remain small. In this per- 



spective, understanding the physical principles of PTR forma- 
tion in multilayers is undoubtedly of importance. Most pre- 
vious works, however, do not really arrive at such principles 
beyond attributing PTR existence to "hidden symmetries" in 
the structures. Instead they draw rather formal conclusions in 
terms of the widely employed transfer matrix method HUO]. 
Such conclusions will benefit from an interpretation to reveal 
their physical meaning. 

In this paper, the question of PTR presence in multilayer 
spectra is addressed from another, more physical than com- 
putational standpoint. Perfect transmission in any multilayer 
(however complex) is seen to be governed by the same princi- 
ples of multiple-beam interference as in a simple Fabry-Perot 
interferometer. Transmission and reflection spectra of any 
multilayer are recovered using recurrent Airy formulas, and 
conditions for any two structures to form PTRs when stacked 
together are derived explicitly. From these conditions, known 
results such as PTRs in mirror-symmetric multilayers natu- 
rally follow. Moreover, it becomes possible to engineer struc- 
tures with PTRs on purpose. As an example, a structure com- 
prising a Fabry-Perot interferometer adjacent to a ID photonic 
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FIG. 1: (Color online) Example transmission spectra of (a) sym- 
metric Cantor multilayer BABAAABAB and (b) non-symmetric Fi- 
bonacci multilayer BABABBAB. A and B correspond to single lay- 
ers with «a = 1 -55, d\ = 76 nm and ;?g = 2.3, d^ = 113 nm, so that 
M A^A = "B^B = 7tc/2a\) = Ao/4 for Xq = 700 nm as in fll5tl . The in- 
sets show an enlarged view of transmission peaks marked by arrows. 



2 



(a) 



(C) 



(b) 

1 

O.i 
<±_ 0.6 
=• 0.4 

0.2 


1.4 

1.2 : 
2 1 
& 0.S 

0.6 



\\ // \\ // \\ « = 15 // 

^1 = 3.5. 


■•V 


■•V 

a? 

A- 


■V 

A' 



1 2 3 4 5 6 

Frequency oj/coq 

FIG. 2: (Color online) (a) A single layer and (b) its transmittance 
| 7a | 2 and phase shift of the reflected wave (p^, given by Eq. 0, for 
different values of «a [^a is chosen in accordance with Eq. 10]; (c) 
an example of a composite S\S2 structure described by Eq. 10. 



crystal is proposed. This structure is shown to feature both 
perfect transmission and highly asymmetric, strongly local- 
ized electric field profile. Localization strength and asymme- 
try can be controlled independently by structure design. 

In SectionHH theoretical background on using Airy-like for- 
mulas for calculating the optical spectra of complex multi- 
layers is given. Section [Til] follows with application of these 
formulas to arrive at the principles encompassing all possible 
cases of PTRs in multilayers. Specific cases such as mirror- 
symmetric and Thue-Morse multilayers are considered, too. 
Section [IV] further employs these principles proposing a de- 
sign for a structure featuring PTRs as well as strongly local- 
ized and highly asymmetric electromagnetic field distribution. 
Finally, SectionMsumm arizes the paper. 



II. RECURRENT AIRY FORMULAS 

We begin by considering a single dielectric layer (labelled 
A) with refractive index n = «a an d thickness d^, located in 
a homogeneous dielectric medium with n — «o (Fig. [2^). Re- 
flection and transmission coefficients of such a layer are given 
by well-known Airy formulas (see, e.g., lfl7l,[l8ll ) 



media labelled by i and j (the wave is incident on the inter- 
face from medium i to medium j). Note that R and T ob- 
tained by Eq. ([TJ are complex and contain information about 
amplitude as well as phase of the reflected and transmitted 
wave. The "usual" intensity-related reflectance and transmit- 
tance are given by \Ra\ 2 an d |7a| 2 , and it can be seen that 
\R-a\ 2 + I 7a 1 2 = L as is obvious from energy conservation. 

In such a simple system, the only frequency-dependent 
quantity is the phase 5a- Since roA = — /"AO and ?oa^ao = 
1 — r A() , it follows that 7\ = 1 whenever 5a = mn for inte- 
ger m. Physically, this corresponds to constructive interfer- 
ence of forward-propagating partial beams inside the layer, 
to occur when its optical thickness is an integer multiple of 
a half-wave. Hence, a single layer features equidistant PTRs 
like a Fabry-Perot interferometer, albeit with poor-quality mir- 
rors (see Fig.lJJ?). The PTR frequencies are 2mcOo with C0q de- 
fined by a well-known quarter-wave (QW) condition 



(n B d B =)n A d A = 7z;c/(2a>o) = Aq/4. 



(2) 



Similarly, let Si and S2 be arbitrary multilayers, e.g., ar- 
bitrary combinations of A and B layers as in Fig. Q] (but of 
course, what follows remains valid way beyond this example). 
Let the reflection and transmission coefficients R$ and T$ be 
known for 5 = 51, S2, and §1 (a bar over Si denotes that Si is 
traversed in the reverse direction). Inserting an infinitely thin 
layer of the ambient medium between the structures (Fig. [2};), 
we can recover the reflection and transmission for the com- 
posite S1S2 multilayer stack: 

p n 1 01 02 Ql t 

K SiS 2 — "Si + "j d — ^—7 ^Si5 2 — 



T Si T s 2 

1 -RsfS: 



(3) 



Note that Eqs. (0 follow from Eqs. (HJ for 5 = 0, and that 
the energy conservation |7j| 2 + \Rs\ 2 — 1 holds. Also note 
that it is critical that both amplitude and phase of R$ and T$ 
is known. In a lossless, linear system one can make use of 
time reversal to relate the spectra of Si and Si as T$ = T$, 
Rg/Ts = -(R S /T s )*. 

By first taking £1,2 to be single layers with reflection and 
transmission spectra given by Eqs. (Q~|i and then using Eqs. (0 
and ([1) in a recurrent fashion, we have a way to calculate 
transmission and reflection spectra for a multilayer of any de- 
gree of complexity. Because such recurrent calculation in- 
volves obtaining reflection and transmission coefficients for 
many intermediate structures, it is numerically less efficient 
than the transfer matrix method. However, the recurrent pro- 
cedure is often adopted for the sake of analytical insight into 
the spectral properties of structures with internal symmetries, 
as was demonstrated, e.g., for fractal multilayers JaS[I3,[3l- 



III. CONDITIONS FOR PERFECT TRANSMISSION 
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(1) 



Here 5a = ( (o/c)n a^a is the phase accumulated by the wave 
in the layer. r i; and fy are Fresnel reflection and transmis- 
sion coefficients, respectively, of an interface between the two 



Our goal is to formulate the existence conditions for a PTR 
in transmission spectrum of an S1S2 stack. From Eq. (0, 
\Ts 1 s 2 (co) \ can be obtained as 



\ T SiS 2 \ = 



\T Sl \\T S2 



1-1**1 1 I 



(4) 



3 



where (p§ and cp$ 2 are the phases of R§ and R$ 2 , respectively. 

Since \R\ 2 + \T\ 2 = 1 in lossless structures, Eq. © can be 
rewtritten in the form 



\T Sl s 2 \ 



1-1*1 



l-\R 2 \' 



|l-|*iN* 



(5) 



<2 



where we have denoted \Ri \ = \R$ l \ = Z?^ , |/??| = |*s 2 |, and 
<p = <p§ + (ps 2 for brevity. If the denominator in Eq. © is 
non-zero, the PTR condition \TsjS 2 \ = 1 is equivalent to 



1-1*1 



1-1*21 



= (I - \Ri\\R 2 \cos<p) 2 + {\Ri\\R 2 \sm<py 



which reduces to 



|*i | 2 + |*2 1 2 = 2|/?i||/? 2 |cosp. 

Obviously, this equation always holds if \R\ \ = |* 2 | = 0, 
which becomes one possible case for PTR, and never holds if 
1*1 1 = 0, |*2 1 or vice versa. In all other cases |*i 1 1* 2 | ^ 
so we obtain 



cos cp = 



l*i 



1*2 1 



2|*l||*2| 



= 1 



l*ll-l*2l) 

2|*l||*2| 



> 1. 



(6) 



If \Ri \ 7^ |/? 2 |, the right-hand side of Eq. (|6]l is strictly 
greater than unity, so no PTR can exist because the condi- 
tion cos cp > 1 cannot be met. If | * 1 1 = |* 2 1 , PTRs can and do 
occur whenever cos cp = 1 . 

For completeness, note that the limiting case when the de- 
nominator in Eq. (0 equals zero results in 

l + (|* 1 ||* 2 |) 2 -2|* 1 ||* 2 |cos(p = 0. 

If \R\ \ \R 2 \ = 0, this equation is false. Otherwise, it can be 
rewritten as 



cos cp = 



1 



l*l| 2 |*2| 2 

2|*i||* 2 | 



= 1 



(1-1*1 1 1*2 1 

2|*i||* 2 | 



>1, 



and can only be satisfied if cos cp = 1 and \Rj 1 1* 2 | = 1. Since 
the reflectance can never exceed unity, the latter implies that 
1*1 1 = l*2| = 1> i- e -> the structure should consist of two per- 
fect mirrors. Such an extreme case causes the right-hand side 
of Eqs. © and © to be indeterminate. This indicates that the 
approach based on the interference of partial waves [Eq. (f3]l] 
becomes invalid with perfect mirrors when there are no par- 
tial waves to interfere. However, this extreme can safely be 
ruled out by assuming that perfect transmission is impossible 
in structures involving perfect mirrors. 

As a result, we have obtained two possibilities for PTR ex- 
istence. First, when \R\ \ = \R 2 \ = 0, or, in the original notation 
ofEq. ©, 



\T Sl 



\Ts 2 



1. 



Second, when \R\ \ = \R 2 \ and cos<p = 1, i.e.. 

\T Sl \ = \T S2 \^L 
<Psi + <Ps 2 = 2m7 *- 



(7) 



(8) 
(9) 



The first condition given by Eq. essentialy means that 
whenever the individual structures Si and S 2 both have a PTR 
at exactly the same frequency, the composite stack S\S 2 will 
always have a PTR at that frequency. In fact, this conclu- 
sion could have been drawn directly from Eq. ©. It is easily 
explained by the fact that if Si and S 2 are both perfectly trans- 
parent, no reflection at the Si /S 2 interface can occur. Hence 
the incident wave is fully transmitted and there is no possibil- 
ity for the reflected wave to form. This is why, e.g., all QW 
multilayers, where all layers conform to Eq. ©, have PTRs at 
CO = 2m(Oo just as any one of the constituent layers. 

The second condition [Eqs. ©-©I is more interesting be- 
'cause it explains how PTRs are formed in the spectral regions 
of the composite structure where there were no PTRs for ei- 
ther Si or S 2 . Indeed, exp[i((pg + <Ps 2 )] — 1 m tne denomi- 
nator in Eq. (0]) renders it equal to the numerator and causes 
\T SlSi \ = 1 although \T Sl \ = \T S , \ ^ 1. The PTR formation 
here can be explained by regarding the composite structure as 
a Fabry-Perot interferometer with very complex mirrors. To 
begin with, the resonance occurs if the interference between 
partial waves is constructive, i.e., if all the partial waves aris- 
ing from multiple reflection are in phase, as given by Eq. (O. 
Then, the resonance is perfect if the mirrors in the interferom- 
eter are balanced and have equal reflectivity [Eq. ©]. 

Eqs. (O-© let us easily see why a mirror-symmetric struc- 
ture readily supports PTRs while most other structures do not. 
Mirror symmetry means = Si, so it is obvious that \Ts t \ — 
\Tg 2 \ and <p^ = (p$ 2 for all frequencies. The only remaining 
condition to be fulfilled is Eq. (O, i.e., <pj = mn. Since 
the phase of the reflected wave varies monotonically between 
transmission resonances in any multilayer with rather few ex- 
ceptions lfl9l Eoll . there should be numerous points where it 
crosses mn (e.g., for one layer it happens for CO — (2m — l)a>o, 
see Fig. |2j>). These points necessarily result in PTRs, as can 
be seen in Fig.(3^. It is seen that for any PTR cos2^ [ = 1, 
except at co = and CO = 2coo where PTRs result from Eq. ^} 
rather than from Eqs. dHJ-©. 

Another simple example would be 52 = Si, i.e., when the 
same structure is repeated twice in the stack. Again we have 
\T$ 2 \ = \T§ \ = 1 7s, I for all frequencies. However, Eq. (O here 
assumes a different form, namely, cps, + % = 2m%, which 
is more difficult to satisfy (compare Figs. [3ji and[3]3). As a 
result, the double-stack structure S\Si exhibits only half as 
many PTRs as does its mirror symmetric counterpart S1S1. 
Both mirror symmetry and stack doubling contribute to PTR 
formation in periodic structures, e.g., one-dimensional pho- 
tonic crystals. Note that if Si is asymmetric, so is S1S1, and 
this case can be regarded as the simplest asymmetric multi- 
layer featuring PTRs. 

Eqs. I©-© also encompass more exotic cases involving 
intrinsically asymmetric structures. Consider Si consisting of 
arbitrarily arranged A and B layers so that «a<^a = »b<^b as 
in Eq. ©, and S2 obtained from Si by substitution A f-> B. 
The resulting structure is very asymmetric (see Fig. [3J;), yet 
it can be shown to feature PTRs. This was observed by Nava 
et al in [J_4] for Fibonacci structures and further pointed out 
by Grigoriev and Biancalana [15] who named such struc- 
tures "Thue-Morse conjugated" because one particular case 
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FIG. 3: (Color online) PTRs in composite structures Si 52: (a) mirror-symmetric (52 = Si); (b) double- stacked (52 = Si); (c) Thue-Morse con- 
jugated. Top: transmission spectra of the constituent structures Si. 2. Middle: transmission spectrum of the whole structure. Bottom: spectral 
dependence of the phase factor cos(<Pj + (p$ 2 ) as in Eq. ([9). The dashed vertical lines show the location of PTRs when both Eqs. <[8j and ((9) 
hold. The dotted lines with arrows designate the peaks that fail to be PTRs due to violation of either Eq. l|9} [in (b)] or Eq. [in (c)]. The 
insets represent a blown-up view of some peaks to determine whether or not they are PTRs. 



of such structures, obtained by repeatedly applying inflation 
rules A — > AB, B — > BA, represents the well-known Thue- 
Morse sequence 112 ill . Fig. [3fc shows calculation results for 
S1S2 with the same Si as for the previous examples (Figs. [3^- 
b). Similarly to these, there are numerous frequencies where 
cos(^ ] + (ps 2 ) = 1 and Eq. (0 is satisfied, and each such 
frequency represents a transmission peak. However, only a 
part of these peaks turn out to be PTRs (see insets in Fig. [3j, 
namely the ones that simultaneously satisfy Eq. ([8]). The rig- 
orous proof of how the fulfillment of these conditions results 
from the Thue-Morse symmetry can be given and is expected 
to appear in a forthcoming publication by Grigoriev et al. 



IV. PERFECT TRANSMISSION 
IN HIGHLY ASYMMETRIC STRUCTURES 

Eqs. I®-® can be employed to engineer a structure of 
any predefined geometry with a PTR at the desired wave- 
length just by varying the refractive index and thickness of 
the layers involved. Indeed, modifying «b/«a in ^2 with- 
out violating Eq. (O changes the value of transmittance and 
reflectance while keeping the phases relatively intact (see 
Fig. [2J3). This aids in fulfilling Eqs. ([8]) and © simulta- 
neously and forms a PTR in the S1S2 structure. Subse- 
quently varying C0q in Eq. (O for both Si and S2 causes all 



the spectra (both amplitude and phase) to scale uniformly, 
thus bringing the PTR to the chosen value of the wavelength. 
Similarly designed dual-interferometer structures of the type 
(AB) m (BA) m (A'B') m (B'A') m were shown to possess PTRs 

Our objective for this paper is to arrive at a design for mul- 
tilayers with highly asymmetric light localization at a PTR so 
as to facilitate the non-reciprocal operation in a nonlinear op- 
tical diode 1 15ll . A straightforward way to achieve the desired 
asymmetry is to stack S\ featuring a strongly localized mode 
with 52 having an extended mode, and to match the frequen- 
cies of the corresponding resonances. 

An obvious choice for Si with a maximally localized mode 
would be a periodic QW multilayer with a half-wave defect, 
or, in other words, a Fabry-Perot interferometer surrounded 
by Bragg mirrors, so that Si = (BA) k (AB) k . If Eq. © holds, 
a sharp transmission resonance occurs exactly at a>o (see 
Fig. |4p). On the contrary, the modes are known to be max- 
imally extended at CO = 2(0q in any QW multilayer. By dou- 
bling the thickness of each layer, this frequency can be halved 
to exactly match the resonance for S\ . A double-periodic pho- 
tonic crystal structure of the type (AABB) m can thus be used 
as ^2- The resulting stack the has the geometry (Fig.|4^) 

S1S2 = (BA) * ( AB ( AABB ) m . (10) 

This design has the obvious advantage that both reso- 
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FIG. 5: (Color online) Electric field localization profile at the PTR 
frequency co = Wq (cf. inset) for the proposed structure design 
[Eq. JLQb l (a) for k = m = 3 as in Fig. |U (b) for k = 3, m = 5 (en- 
hanced asymmetry), and (c) for k = 5, m = 3 (enhanced localization 
strength). 



FIG. 4: (Color online) (a) The proposed design of a highly asymmet- 
ric multilayer featuring PTRs given by Eq. l |10t , along with transmis- 
sion spectra for (b) S\, (c) 52, and (d) S\S%, 



nances in question are PTRs (and they are exactly frequency 
matched), so there is no need to go as far as Eqs. ©-(HI) and 
the resulting PTR in S\Si is assured due to Eq. The sec- 
ond advantage is that frequency matching always occurs at 
A = Ao/4 so it is easy to design the structure for any desired 
wavelength using any materials at hand. Fig.HJl confirms the 
existence of a PTR, as does explicit numerical calculation of 
|Tsift(<Wo)|> yielding 1 within limits of machine accuracy. 

The electric field intensity distribution at the resonant fre- 
quency is shown in Fig. [5] It is seen that the localization is 
clearly asymmetric and mainly present in Si. Using the same 
materials as in Ref. 11511 , comparable localization strength is 
observed for a structure about 3 times thinner and having 32 
layers instead of 64 (see Fig. [5}:). The design of Eq. ( fTOb also 
allows to control the localization strength (by varying k) and 
the asymmetry (by varying m) independently and in a wide 
range for a relatively minor change to the number of layers 
(compare Fig. [5^— c). This is opposed to changing the num- 
ber of generations in a Thue-Morse sequence, which would 
double or halve the number of layers at once. A possibility to 
build PTR-enabled structures with desired localization prop- 
erties using relatively few layers is important from the prac- 
tical point of view because losses would obviously be more 
detrimental to perfect transmission in thicker structures IU4I1 . 

Note, finally, that the choice of geometry for S2 is rather 
arbitrary as any arrangement of AA and BB will produce the 
same extended-mode PTR at coq. This choice of geometry can 
be regarded as an additional design tool to influence the trans- 
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FIG. 6: (Color online) Transmission spectrum of the structure in 
Fig.f4^ for ;ia = 1.55, jib = 2.3 as in flal (dotted) and for increased 
n B /"A by setting = 1 (solid). 



mission spectrum around the PTR. For a periodic geometry 
52 = (AABB)'" used in Eq. ( fTUI ). the two band gaps around 
(1 ± 1/2) (»o brought about by 62 (see Fig. |4};) can overlap 
with the gap around ©0 for Si (Fig. Eh). This would widen 
the region of predominantly low transmission surrounding the 
designed PTR, which can prove useful. For the materials 
adopted throughout this paper from Ref. 11511 . it is not yet the 
case, but the gap overlap can be achieved by increasing «b /«a 
(Fig. [6]). It is seen that the PTR then becomes very isolated in 
the transmission spectrum. 
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V. CONCLUSIONS 

Using the formalism of recurrent Airy formulas, the condi- 
tions for a mulilayer structure to exhibit perfect transmission 
resonances [Eqs. (0-(|9]l] are formulated rigorously in such a 
way that possibilities for PTRs can be directly envisioned at 
the stage of structure design. Following the previous results 
lfl4ll . it was shown that mirror symmetry is a sufficient but not 
necessary condition for PTR existence. PTRs are shown to 
be possible in asymmetric structures, including Thue-Morse 
conjugated multilayers II 1 511 . Based on frequency-matched 
PTRs in structure parts, the design for a combined Fabry- 
Perot/double-period photonic crystal multilayer was proposed 
(Fig. |4^). This structure was shown to feature perfect trans- 
mission resonances with strongly localized and highly asym- 
metric spatial distribution of electric field intensity (Fig. [5]). 
Strength and asymmetry of localization can be controlled in- 
dependently by changing the design parameters, keeping the 



number of layers reasonably small. It is expected that mul- 
tilayers of this kind would enhance non-reciprocal transmis- 
sion if they contain nonlinear materials, improving the perfor- 
mance of optical diodes and similar devices. 
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